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I.  Research  Sunmary 

Our  research  goal  has  been  to  investigate  the  integration  of  the  Calspan 
(2),  Bowman  (3),  and  Linkman  or  Bubble-man  (4)  programs  to  provide  a visual 
(graphic)  facility  for  examining  and  evaluating  these  models.  This  research 
included  transforming  the  Calspan  body  position  data  to  the  Linkman  format, 
combining  the  body  model  with  a given  cockpit  model,  displaying  both  model  and 
cockpit  together  in  a readily  visualizable  manner,  and  detecting  collisions 
between  the  body  and  the  cockoit.  Each  of  these  areas  is  sunmarized  below, 
further  details  are  provided  in  the  Appendices. 

The  Calspan-to-Bubble-nvan  conversion  program  interfaces  these  two  systems 
by  converting  data  frcm  one  of  the  Calspan  output  files  into  a suitable  input 
file  for  the  Bubble-man  program.  By  using  files,  no  direct  interaction  or 
modification  of  either  program  is  required.  The  significant  portion  of  this 
effort  is  based  on  locating  the  appropriate  data  in  the  Calspan  file  and  con- 
verting coordinate  data  from  one  model  into  the  (different)  model  in  Bubble-man. 
We  used  the  body  segment  angular  displacements  to  compute  the  joint  angles 
required  in  Bubble-man. (Details  are  given  in  Appendix  1.)  Figure  1 shows 
four  positions  of  the  body  model  which  were  produced  by  the  Calspan  simulator 
and  displayed  with  the  Bubble-man  program. 

Although  the  Calspan  program  is  capable  of  writing  out  its  variable  in- 
formation at  as  small  a time  step  as  desired  (but  of  course  not  be lew  the 
integration  step  size),  it  still  may  be  desirable  to  interpolate  between  the 
Calspan  results.  Interpolation  might  be  needed  for  producing  an  animated  film 
sequence,  for  instance.  For  such  uses,  an  APL  program  was  developed.  This 
program  reads  a Bubble-nan  angle  file,  such  as  that  produced  by  the  Calspan- 
Linkman  interface  program,  and  creates  a new  file  of  angles  at  a shorter  time 
interval.  Hew  much  shorter  is  determined  by  an  input  parameter.  The  inter- 
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polation  is  a straight- forvard  linear  interpolation.  The  APL  listing  is  given 
in  Appendix  2. 

To  produce  understandable  pictures  of  the  human  body  model  and  cockpit 
data  from  the  Bcwman  program,  techniques  were  developed  for  combining  Bubble-man 
with  planar  polygon  cockpit  surfaces  (Figure  2).  Since  Bubble-man  consists  of 
overlapping  spheres,  no  existing  techniques  could  be  used  directly.  Instead, 
the  cockpit  polygons  and  the  Bubble-man  spheres  are  entered  into  two  arrays 
of  points:  one  representing  the  actual  intensity  of  that  point  on  a raster 
display;  and  the  other  representing  the  depth  (distance  from  the  observer)  of 
the  closest  scene  component  at  that  point.  As  each  Bubble-man  sphere  or  planar 
polygon  is  processed,  it  is  checked  against  the  current  depth  values  at  each 
point  it  affects.  If  its  depth  is  greater,  then  nothing  is  changed.  Otherwise 
the  depth  value  is  changed  to  the  closer  point  and  the  intensity  array  is  up- 
dated as  well.  The  result  is  a shaded,  solid  rendering  of  body  and  cockpit 
(Figure  3). 

Given  the  body  position  and  cockpit  model,  collisions  (intersections)  be- 
tween them  may  be  discovered  by  simply  computing  the  intersection  between  each 
sphere  and  polygon.  In  practice  (since  there  are  many  spheres  and  polygons) 
the  number  of  comparisons  is  reduced  by  checking  each  body  segment  against  each 
polygon,  then  checking  each  sphere  in  the  segment  if  an  intersection  were 
possible  by  the  first  test.  Further  details  and  an  example  of  the  collision 
detection  process  are  given  in  Appendix  3. 

Further  information  on  the  Bubble-man  model  and  its  properties  may  be 
found  in  Reference  (1).  We  believe  that  the  proposed  integration  of  the 
modelling  and  simulation  systems  can  be  readily  achieved  as  demonstrated  by 
our  results,  and  that  such  an  integration  is  desirable  in  order  to  better 
utilize  and  visualize  human  biodynamic  information. 


Figure  2a.  Cockpit  drawn  with  lines  on  raster  display. 


at  drawn  with  hidden  surfaces  removed. 


Figure  3.  Mi  min  body  model  in  cockpit. 
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II.  Technical  Reports 

Reference  (1)  above  shall  serve  as  a manuscript  description  of  the  overall 
thrust  of  this  project.  Other  relevant  details  are  provided  as  Appendices  to 
this  final  report. 

Ill .  Publications 

Reference  (1)  has  been  submitted  for  publication.*  The  cognizant  Program 
Director  for  the  Office  of  Naval  Research  will  be  notified  and  supplied  with 
reprints  upon  acceptance. 

IV.  Conclusions 

We  have  demonstrated  that  a solid  model  of  the  human  body  (Bubble-man) 
may  be  positioned  by  joint  and  segment  data  obtained  from  an  existing  movement 
simulator  without  modification  of  the  latter.  We  have  further  demonstrated 
that  the  body,  once  positioned,  may  be  placed  in  a planar  polygon  environment 
such  as  used  for  cockpit  design  and  evaluation.  Besides  efficiently  creating 
a graphics  display  of  the  composite  (body  plus  cockpit),  collisions  between 
the  body  and  itself  or  the  cockpit  may  be  easily  computed,  yielding  concise 
(textual)  descriptions  of  contact  or  intersection  points.  Ibis  information  is 
essential  to  the  evaluation  of  the  cockpit  design  on  the  human  occupant. 

By  integrating  these  three  systems  we  open  the  possibility  for  video 
recording  and  preservation  of  real  or  simulation  data:  movements  of  the  body 
may  be  examined  as  they  occur  in  the  cockpit  or  vehicle.  Computer  graphics 
permits  the  placement  of  the  "observer"  at  any  position,  thus  multiple  views 
of  the  same  movement  sequence  are  easily  produced.  The  full  three-dimensionality 
and  solidity  of  the  model  also  provide  a degree  of  realism  unequalled  in  other 
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models.  Such  animations  could  be  used  to  provide  long-term  archives  of  move- 
ment data,  independent  of  computers  or  programs.  By  suitable  video  mixing, 
real  and  simulated  data  could  be  overlaid  for  validation  of  the  simulator  itself. 
V.  Major  Accompli shment s 

We  may  itemize  our  conclusions  as  follows: 

1.  Bubble-man  may  be  positioned  by  data  obtained  from  existing  movement 
simulators. 

2.  The  Bubble-man  model  effectively  portrays  a body  position  on  a raster 
graphics  display:  hidden  parts  are  removed;  the  depth  relationships 
are  apparent;  and  joint  deformations  are  eliminated. 

3.  Acceptable  graphic  images  of  Bubble-man  may  be  produced  on  vector 
drawing  as  well  as  raster  displays. 

4.  Joint  positions  obtained  from  a simulator  may  be  interpolated  to 
achieve  additional  (intermediate)  positions  for  an  animated  display. 

5.  A movement  simulator  may  be  used  to  generate  body  positions,  and  that 
body  may  then  be  placed  in  a planar  polygon  environment  (such  as  an 
airplane  cockpit).  Both  may  be  displayed  with  hidden  parts  removed 
and  visible  parts  shaded. 

6.  Collision  detection  between  Bubble-man  and  a planar  polygon  environment 
is  simple  and  efficient,  yielding  concise  (textual)  descriptions  of 
contact  or  intersection  points. 

7.  Demonstration  of  the  feasibility  of  transferring  simulation  data  to 
video  tapes  of  the  body  model  in  a given  environment,  useful  for 
evaluation  and  archival  storage  of  experimental  data. 
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Appendix  1 

Introduction 

This  appendix  is  concerned  with  the  details  of  interfacing  three  different 
computer  programs:  The  Calspan  Crash  Victim  Simulator,  The  Bowman  Cockpit 
Geometry  Program,  and  Linkman  or  Bubble-nan.  Each  program  is  useful  in  its 
own  right,  but  used  together  they  can  become  powerful  tools  for  the  analysis 
and  validation  of  biodynamic  models. 


Calspan  to  Linkman  Conversion 

The  main  purpose  of  the  Calspan  to  Linkman  Conversion  program  is  to  act 
as  an  interface  between  two  existing  programs:  The  Calspan  Crash  Victim 
Simulation  and  The  Linkman  or  Bubble-nan  program.  The  interface  program  (called 
CALBUB)  accomplishes  this  by  reading  data  from  one  of  the  Calspan  output 
files  and  preparing  an  input  file  for  the  Linkman  program.  It  works  entirely 
from  files:  There  is  no  direct  interaction  with  either  the  Calspan  or  Linkman 
programs. 

It  may  seem  at  first  sight  that  the  interfacing  task  is  trivial,  since 
the  Linkman  needs  joint  angles  as  input,  and  the  Calspan  program  outputs 
joint  angles  (if  requested).  However,  the  joint  angles  output  by  the  Calspan 
program  are  intended  to  explicate  the  joint  torques,  not  describe  the  position 
of  the  body.  Therefore,  they  are  given  as  flexural  and  twist  angles  from  the 
"principle  axes"  defined  for  the  joints.  It  is  not  possible  to  convert  these 
flexural  and  twist  angles  into  the  yaw,  pitch,  and  roll  angles  needed  by  the 
Linkman  program.  This  is  because  the  flexural  and  twist  angles  only  represent 
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two  degrees  of  freedom,  the  yaw,  pitch,  and  roll  system  expresses  the  full 
three  degrees  of  freedom. 

In  contrast  to  these  difficulties,  the  segment  angular  displacements 

! 

represent  a reliable  source  of  data  for  the  joint  angles.  Normally  the  segment 
angular  displacements  for  all  segments  are  written  out  by  the  Cal span  program, 
and  the  conversion  into  joint  angles  is  conceptually  clean,  but  by  no  means 
trivial.  The  program  CALBUB  uses  these  segment  angular  displacements. 

The  program  will  be  described  in  three  parts:  the  reading  of  the  Calspan 
output  file,  the  angle  conversion,  and  the  writing  of  the  Linkman  input  file. 
Reading  of  Calspan  File 

The  Calspan  Simulator  creates  a number  of  output  files.  One  of  them,  TAPE4, is 
designed  for  time-curve  plotting,  and  contains  the  angular  displacements  of 
the  segments  (if  requested).  In  order  to  read  this  file  and  extract  only  the 
segment  angular  displacements,  it  is  necessary  to  understand  the  precise  file 
structure  f nd  positions  of  the  data. 

!,  TAPE4  consists  of  a number  of  records,  each  record  of  which  consists  of 

all  the  variables  at  a particular  time.  Thus  the  overall  structure  of  the  file 

i is  as  in  Figure  1. 

I , , , , 


all  variables  at  tg 

all  variables  at  t^ 

all  variables  at  t 
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Within  each  time  record,  the  variables  are  divided  into  ten  groups,  named  as 


follows: 

LA 

linear  acceleration 

LV 

linear  velocity 

LD 

linear  displacement 

AA 

angular  acceleration 

AV 

angular  velocity 

AD 

angular  displacement 

JP 

joint  parameters 

PS 

plane-segment  contact  faces 

BS 

belt-segment  contact  faces 

SS 

segment- segment 

We  will  call  each  of  these  groups  a code,  and  will  refer  to  them  by  the  above 
two  letter  designations.  For  each  code,  there  may  be  a number  of  different 
data  sets.  These  data  sets  can  usually  be  identified  with  segments,  but  not 
always.  For  example,  for  the  code  LA,  there  may  be  three  data  sets:  neck 
c.g.,  head  c.g. , and  head  at  neck  pivot.  In  this  example,  there  are  two  data 
sets  for  the  same  head  segment.  The  last  level  of  structure  is  represented  by 
components . Each  data  set  has  a number  of  components,  for  example,  x,y,z,  and 
resultant.  The  number  of  components  is  determined  by  the  code,  not  the  data 
set.  Thus  all  LA  data  sets  always  have  4 components  each,  all  JP  data  sets 
have  6 components,  etc. 

To  summarize,  the  structure  of  the  Calspan  output  file  can  be  characterized 
by  4 different  levels: 

1.  records 

2.  codes 

3 . data  sets 

4 . components 

In  order  for  the  interface  program  CALBUB  to  read  the  Calspan  file,  it 
must  be  told  the  particular  file  structure  to  expect.  For  this  purpose,  the 
CALBUB  program  reads  file  definition  information  from  TAPE5.  This  file 
(TAPE5)  is  the  only  input  to  CALBUB  other  than  TAPE4,  the  Calspan  output  file. 
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We  have  now  set  enough  background  to  discuss  the  precise  form  of  TAPES,  the 
input  definition  file. 

CARD1.  either  DEBUG  or  NO  DEBUG  (format  A10) 

This  controls  the  printing  of  debugging  information  in  TAPE2. 

CARD?.  The  number  of  data  sets  per  code  (free  format) 


CARD3.1 

; 

CARD  3.  NANGDS 


The  names  of  the  data  sets  (i.e.  segments)  for  angular 
displacement  data  (format  A10) 


CARD4.  number  of  joints 


CARD5.1 

0 

CARD5.JMAX 


The  names  of  each  joint  and  its  two  associated  segments 
(format  3A10) 


Card 2 specifies  the  file  structure  by  listing  the  number  of  data  sets 
per  code.  Cards  3 give  the  names  of  the  data  sets  for  the  angular  displacement 
data  which  we  are  interested  in,  and  Cards  5 associate  these  segment  or  data 
set  names  with  the  joint  names.  A Sample  14  segment  input  deck  is  shown  below. 


NO  DEBUG 

:L  1 y 1 » :i.  2 y : 

i.  r 

GT 

U I 

N 

H 

ELIL 

! r 2 * 1 

:i  y 4 r 0 f 0 * 0 

El.  t. 

EE 

L III.. 

LI  l 

L 1 

ELIA 

EL.A 

UJA 

L.L.A 

1.4 

P 

LT 

GT 

+ 

U 

GT 

LIT 

* 

NP 

LIT 

N 

* 

HP 

N 

H 

* 

EH 

LT 

ELIL. 

* 

EK 

ELIL 

ELL 

* 

EA 

ELL 

EF 

* 

LH 

LT 

L.LIL 

* 

LK 

LLJL 

Ll  L. 

* 

LA 

LL.L 

LF 

* 

EG 

LIT 

ELIA 

* 

EE 

ELIA 

E!  L 

* 

L.S 

LIT 

L LI  A 

* 

L.E 

L..UA 

L 1 . 1 

# 

T 

X... 
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After  the  input  file  TAPE5  is  read  in  by  subroutine  RDEF,  a routine  PREPRT 
is  called  to  prepare  for  the  reading  of  the  Calspan  tape.  This  routine  uses 
the  definition  information  to  compute  a mapping  vector  POINTER  which  specifies 
exactly  which  positions  within  each  time  record  of  the  Calspan  file  should  be 
stored  for  future  use.  Then  subroutine  RTAPE  is  called,  and  this  reads  in 
each  record  and  stores  away  the  appropriate  information  according  to  the 
mapping  vector.  When  RTAPE  is  finished,  all  the  Calspan  angles  are  stored 
within  the  array  CALANG  as  follows: 

‘o  »]  P]  rl  y2p2r2  ••••ynpnrn  *1  yl  P1  rl  "" 

The  angles  are  stored  in  the  order  they  appear  on  the  original  Calspan  tape, 

i.e.,  the  order  represented  by  cards  3 of  the  input  definition  file. 

This  represents  the  completion  of  the  first  phase  of  the  CALBUB  interface 

program:  the  Calspan  tape  has  been  read  and  the  segment  angular  displacement 

data  stored  in  the  array  CALANG. 

Angle  Conversion 

Before  we  can  discuss  the  conversion  between  the  Calspan  and  Linkman  angles, 
it  is  necessary  to  rake  clear  the  definition  of  the  angles.  Both  the  Calspan 
and  Linkman  coordinate  systems  are  all  right  handed,  with  the  three  angles 
defined  as: 

yaw:  + x to  y 

pitch:  + z to  x 

roll:  + y to  z 

However,  all  the  Calspan  coordinate  systems  have  their  z-axis  pointing  downwards, 
(with  gravity)  whereas  all  the  Linkman  coordinate  systems  have  their  z-axis 
pointing  upwards  (against  gravity).  The  Calspan  convention  is  responsible  for 
an  illusion  that  the  coordinate  systems  in  the  upper  and  lower  body  parts  are 
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different.  They  are  not.  It  is  just  that  the  segments  sometimes  "hang"  under 
the  coordinate  system,  and  sometimes  sit  on  top  of  the  coordinate  system. 

In  order  to  change  the  angles  relative  to  the  Calspan  inverted  coordinate 

system  into  angles  relative  to  an  upright  coordinate  system,  the  following 

/ 

transformatiph  is  performed: 

1.  change  sign  of  yaw 

2.  change  sign  of  pitch 

3.  leave  roll  unchanged 

Once  this  transformation  is  performed,  we  no  longer  need  worry  about  differences 
between  coordinate  system  definitions. 

All  the  Calspan  segment  angular  displacements  are  expressed  as  rotations 
from  the  vehicle  to  the  segment.  The  angles  needed  for  the  Linkman  program 
are  the  joint  angles:  rotation  from  one  segment  to  the  adjacent  segment  con- 
nected by  the  joint.  Hie  first  segment  is  more  proximal  to  the  body  center, 
the  second  segment  more  distal.  The  two  segments  which  together  define  a joint 
angle  are  specified  by  the  card  group  5 in  the  input  definition  file. 

The  method  of  computing  the  Linkman  joint  angles  given  the  two  associated 
Calspan  segment  angles  (expressed  in  upright  coordinate  systems)  is  as  follows. 
From  the  segment  1 angles,  a rotation  matrix  is  computed  which  transforms  from 
segment  1 to  the  vehicle.  Likewise  a matrix  is  found  which  transforms  from 
the  vehicle  to  segment  2.  The  composition  of  these  two  rotations  gives  a 
transformation  from  segment  1 to  segment  2 (See  Figure  2). 


-i*. 


segment  1 


segment  2 
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Figure  2 


The  rotation  matrices  are  direction  cosine  matrices , and  are  computed  as: 


cos  p cos  y 

cos  p sin  y 

-sin  p 

-sin  y cos  r 

cos  y cos  r 

cos  p sin  r 

+sin  r sin  p cos  y 

+sin  r sin  p sin  y 

+sin  y sin  r 

-sin  r cos  y 

cos  p cos  r 

+cos  r sin  p cos  y 

+cos  r sin  p sin  y 

combined  transformation  M^' 

T 

has  this  form.  Once  we 

have  completed 

this  matrix,  our  next  task  is  to  find  three  angles  which  could  give  rise  to 
the  matrix.  These  angles  are  the  Linkman  joint  angles  we  are  seeking.  By 
inspection  of  the  matrix  elements  we  can  derive  the  following  formulas: 
yaw  = arctan  (M(l,2)  / M(l,l)) 
pitch  = arcsin  (M(l , 3 ) ) 
roll  = arctan  (M(2,3)  / M(3,3)) 


-18- 

By  the  nature  of  the  inverse  trigonometric  functions,  these  formulas  leave 
some  ambiguity.  In  particular,  we  could  also  have: 
yaw  + 180 
supplement  (pitch) 
roll  + 180 

This  gives  us  eight  possible  solutions. 

Although  there  may  be  a method  of  choosing  among  these  8 possibilities 
analytically,  we  have  found  it  easier  to  simply  substitute  each  angle  set  back 
into  the  definition  of  the  direction  cosine  matrix  and  see  if  indeed  the  angles 
yield  the  original  matrix.  The  result  of  this  substitution  is  that  two  angle 
sets  are  always  found  which  will  give  the  original  matrix.  These  two  angle 
sets  both  represent  the  exact  same  angular  orientation.  The  set  which  is  closer 
to  that  used  as  the  solution  for  the  previous  time  step  is  chosen  as  the 
solution  for  the  current  time  step. 

At  the  completion  of  the  angle  conversion  phase  of  the  program,  all  the 
joint  angles  have  been  computed  by  the  above  process  and  stored  into  an  array 
BUBANG,  whose  structure  is  exactly  similar  to  the  structure  of  CALANG. 

Writing  of  Angle  File 

The  last  phase  of  the  CALBUB  program  is  to  write  the  joint  angles  on  a 
file  in  a form  ready  to  be  read  in  by  the  Linkman  program.  This  is  rather 
straightforward,  but  there  are  two  further  changes  made  to  the  data. 

First,  the  Calspan  and  Linkman  programs  differ  in  their  definition  of 
ankle  angles.  In  the  Linkman  program,  all  zeros  for  the  ankle  angles  gives 
the  normal  relaxed  standing  position;  in  the  Calspan  program,  all  zeros  means 
the  feet  are  pointing  straight  downwards.  Thus  90  degrees  are  added  to  the 
pitch  of  the  Calspan  ankle  angles  to  be  consistent  with  the  Linkman  conventions. 
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Second,  it  Is  normally  the  case  that  the  Calspan  output  does  not  give 
sufficient  information  for  determining  every  Linkman  joint  angle.  The  simula- 
tion may  not  lave  used  as  many  segments  as  there  are  in  the  Linkman  model,  or 
not  all  of  the  information  may  lave  been  directed  to  the  Calspan  output  file. 

A list  of  all  the  joint  names  'used  in  Linkman  is  therefore  compared  against 
the  joints  define  by  the  Calspan  data  (specified  by  card  group5  on  the  input 
definition  file),  end  for  any  joints  which  are  not  defined,  default  values  are 
substituted.  These  default  values  are  set  in  a data  statement  within  subroutine 
WF1LE. 

With  these  two  changes  made,  the  joint  angles  are  written  out  on  a file 
(TAPF6)  which  can  serve  as  the  input  file  (TAPE4)  for  the  Linkman  program. 

Using  the  Interface  Program 

There  are  three  steps  in  any  use  of  the  CALBUB  program: 

1.  Run  the  Calspan  Simulator  and  save  TAPE4. 

2.  Run  CALBUB  using  job  control  cards  equivalent  to  the  following: 

GET  (TAPE4  = [file  name  of  Calspan  output  file]) 

GET  (TAPE5  = [file  name  of  input  definition  file]) 

GET  (LGO  + CALBUBB) 

FILE  (TAPE4,  RT=S. . . ) 

LDSET(PRESET=  ZERO,  MAP=S,  FILES=TAPE4) 


4 


1 


LGO. 

REPLACE  (TAPE6  = [file  name  for  Linkman  angles]) 
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3.  Run  the  Linkman  program  using  job  control  cards  equivalent  to: 
GCT  (TAPL3  = [file  name  of  man  data]) 

GET  (TAPE4  = [file  name  of  angles]) 

GCT  (PL10UB/  UN= SYSTEM) 

GCT  (UJKMANB) 

LDSET  (PRESCT=ZERO,  MAP=S,  LIB=PL10LIB) 

LNKMANB. 
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Appendix  2 

APL  angle  interpolation  program. 
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Sphere  - polygon  collision  algorithm. 

1.  Find  plane  containing  polygon. 

2.  Find  distance  between  sphere  surface  and  polygon  plane.  If  this 
distance  is  greater  than  zero,  sphere  and  plane  do  not  intersect; 
otherwise  sphere  intersects  plane  in  a (possibly  degenerate)  circle. 

3.  Coordinatize  plane  of  polygon , ‘lien  check  if  circle  intersects  or  con- 
tains enclosing  box  of  polygon.  If  not,  sphere  and  polygon  do  not 
intersect. 

4. a If  center  of  circle  is  inside  polygon,  then  sphere  and  polygon  intersect, 
b.  If  center  is  outside  polygon  and  any  edge  of  polygon  intersects 
circle,  then  sphere  intersects  polygon. 

In  case  the  circle  in  4a  is  a single  point  the  intersection  may  be  called 
a "touch"  (tangency).  Otherwise  the  volumes  overlap.  Also,  in  case  4b,  one 
rrwy  wish  to  distinguish  the  case  in  which  the  circle  is  tangent  to  the  polygon, 
but  not  otherwise  intersecting  it. 

An  example  of  this  algorithm  on  the  body  and  cockpit  shown  in  Tigure  2 
of  the  report  yields  the  following  collision  description.  (The  polygon  names 
are  supplied  from  the  cockpit  data  file,  the  body  segment  names  from  the 
Bubble-man  program. ) 
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